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Collisionless Magnetic Reconnection in a Five-Moment 
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Abstract. We simulate magnetic reconnection in electron-positron (pair) plasma 
using a collisionless two-fluid model with isotropic pressure. In this model the 
resistive, Hall, and electrokinetic pressure terms are absent from the curl of 
Ohm's law, leaving the inertial term alone to provide for magnetic reconnec- 
tion. Our simulations suggest that for pair plasma simulated with isotropic 
pressure fast reconnection does not occur without the aid of sufficient (numer- 
ical) diffusion. We contrast this result with simulations and published results 
showing fast reconnection for collisionless two-fluid plasma with isotropic pres- 
sures and non-canceling mass-to-charge ratios, where Hall effects are present 
and numerical diffusion is small, and with published PIC studies of pair plasma 
which observe fast reconnection and attribute it to nonisotropic pressure. 



1. Introduction 

An important issue of controversy in the magnetic reconnection community 
is the minimal conditions required for fast magnetic reconnection to occur in a 
plasma and the minimal modeling requirements to resolve it [7] . The first attempt 
to model reconnection was carried out by Sweet [9] and Parker [6], who used a 
resistive magnetohydrodynamic (MHD) model. Their approach was only successful 
in modeling a slow form of reconnection and not the much faster reconnection that 
is observed in laboratory and space plasma. 

The Geospace Environmental Modeling (GEM) Reconnection Challenge prob- 
lem was introduced in [3] and studied with a variety of models to identify the es- 
sential physics required to model collisionless magnetic reconnection. This original 
GEM article concluded that all models that include the Hall term in the general- 
ized Ohm's law produced essentially indistinguishable rates of reconnection. The 
only other model in their study that admitted fast reconnection was MHD with 
large anomalous (e.g., current-dependent) resistivity, although as expected it did 
not exhibit the quadrupole out-of-plane magnetic field pattern that appears to 
characterize models which incorporate the Hall term. 
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Since inclusion of Hall effects had been identified as the critical ingredient to 
admit fast reconnection, Bessho and Bhattacharjee studied electron-positron plas- 
mas, for which the Hall term is zero|ll[2]- Their particlc-in-cell (PIC) simulations 
of the GEM problem (with the mass ratio reset from 25 to 1) exhibited fast recon- 
nection for temperature ratios of 5 and 1. They found that the structure of the 
reconnection region still shows an X-point, although as expected the out-of-plane 
magnetic field does not show the quadrupolc structure that appears in models that 
incorporate Hall term effects. 

This prompted us to ask if it is possible to get fast reconnection if the Hall term 
is absent and the pressures are modeled as isotropic. Therefore we chose to study 
reconnection in an electron-positron plasma using a two-fluid model with isotropic 
pressure. 

2. Physical models 

2.1. Particle-in-cell description. Plasma consists of charged particles in- 
teracting with the electromagnetic field. In the absence of gravity and quantum- 
mechanical effects the particles of a plasma satisfy Maxwell's equations and the 
Lorentz force to govern particle motion: 

dtB = -V X E, V • B = 0, 

dtE ^c^\/ X B - J/eo, V • E = a/eo, 

dt-Vp = ^ (|e(xp) + Vp X B(xp)^ , dtXp = Vp, 

J = ^S'p(xp)(7pVp, (J ^^Sp{xp)qp; 

p p 
here B is magnetic field, E is electric field, c is the speed of light, eq is electric 
permittivity, p is particle index, Xp(t) is particle position, Vp(t) = 7pVp is (proper) 

— 1/2 

particle velocity, where 7 = (l — (ij/c)^) w 1 is the Lorentz factor, qp is particle 
charge, rup is particle mass, a is charge density, J is current density, and S'p(x — Xp) 
is particle charge distribution (e.g., a unit impulse function). 

Particle-in-cell (PIC) codes model plasma by attempting to evolve in a mixed 
Eulcrian-Lagrangian framework: particles are treated in a Lagrangian way, while 
the electromagnetic field sits on an Eulerian computational mesh. 

2.2. Boltzmann description. The Boltzmann model replaces the particles 
with particle (probability) density functions /s(x,v,t) of space, velocity, and time, 
for each species s. The Boltzmann equation asserts conservation (or balance) of 
particles in phase space: 

dtfs + Vx • (v/,) + -Vv • (^(E + V X B)fs) = 

here v = 7V « v is (proper) velocity in phase space, and is a collision operator 
which is a function of {v 1— *■ /p(t, x, v)}p, where p ranges over all species. The colli- 
sionless Boltzmann equation (alias Vlasov equation) asserts that Cs = 0. Maxwell's 
equations are coupled to the Boltzmann equation by the relations 
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2.3. Two-fluid model. Multiplying the Boltzmann equation by powers of 
velocity and integrating over velocity space yields fluid equations. Generic two- 
fluid equations for a two-species plasma are: 

(2.1) 
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The variables are defined as follows: i and e are ion and electron species indices; 
for species s € {i,e}, Qs = ±e is particle charge, TOs is particle mass, is particle 
number density, ps = rUsUs is mass density, cts = Isns is charge density, 3s = UgCTs 
is current density, Ps is the pressure tensor, £s is gas-dynamic energy, qs is the 
heat flux, R.^ = —Re denotes the interspecies drag force on the ions, Qr^s denotes 
heating due to friction (drag), and Qi — —Qe denotes the interspecies thermal heat 
transfer to the ions. 

2.4. Collisionless isotropic closure. To close the system we must posit 
constitutive relations for the nonevolvcd quantities. In a collisionless model we 
neglect the terms that come from the collision operator: q^, R^, Qr^s, and Qs- In 
an isotropic model we assume that the pressure tensor is a scalar pressure times 
the identity tensor: — p^I; this leads to the constitutive relation Eg = {3/2)ps + 
Psull2. 

2.5. Ohm's law. Multiplying the momentum equations of each species by its 
charge to mass ratio and summing gives a balance law for net current. Invoking the 
assumption of quasincutrality (cr « 0) and solving this law for electric field gives 
the generalized Ohm's law, 

E = B X u + E', 

where u is the mass-averaged fluid velocity and where the electric field in the frame 
of reference of the fluid is the sum of four terms: 



E'= 77 -J 



-J X B 



iv • (meP. - m,Pe) 

p 

fhifhe 



p 



ft:J 



uJ + Ju 



p 



(resistance) 
(Hall term) 

(pressure term) 

-Jj) ) (inertial term). 



^ We remark that our collisionless model with isotropic pressure seems not to correspond to 
any general physical regime of plasma (although it may apply to particular configurations). A 
model is considered to be physical if it agrees with a physical regime in some physical limit. We 
assume an isotropic pressure tensor yet no resistivity, but for an electron-positron plasma the time 
scale over which particles thermalize is the same as the time scale over which resistive drag force 
seeks to equilibrate the velocities of the two species. 
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Here we adopt the convenction that m,, := rrisje^ and we have assumed that the 
electrical resistance equals 77 • J (where 77 is the resistivity), a simple 

function of the drift velocity, i.e., of the current. 

Ideal MHD assumes that E' = 0, and resistive MHD assumes that E' = 77 • J. 

Substituting Ohm's law into Faraday's law dtlB + V x E = yields 9tB + V • 
(uB + Bu) = V X E', which implies that the flux of B through a surface convected 
by u can only change if the curl of E' is nonzero. 

In an electron-positron plasma the masses of ions and electrons are identical 
and the Hall term vanishes. Also, if pressure is isotropic and density varies slowly, 
then the curl of the pressure term is zero. In the collisionlcss two-fluid model the 
resistivity is zero. So for our two-fluid model reconnection could only happen by 
means of the inertial term (or numerical diffusion). 



3. The five-moment two-fluid model 

The collisionlcss two-fluid equations we solved were 
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These equations were studied extensively by Shumlak and Loverich 2] and Hakim, 
Shumlak, and Loverich [5]. A version of this model with anisotropic pressure was 
also considered by Hakim [4]. 

This system is identical in apearance with the two-fluid system (j2.ip - (|2.2p if 
the correction potentials ip ^^nd (f), which we have added for numerical divergence 
cleaning purposes, are zero. These equations imply a wave equation that propagates 
the divergence constraint error at the speed cxu We select x = 1-05. 

We nondimensionalized this system by choosing typical values of magnetic field 
Bq, (ion) number density tiq, particle charge qo = e, and combined particle mass 
mo = m.i + rrie- This implies a choice (1) of characteristic time scale w"^ := 
where Ug is the cyclotron frequency of a typical particle, (2) of characteristic velocity 
VA ^^amona ' ^ typical Alfveu speed, where /xq := (c^eo)^^ is the permeability 
of magnetic field and (3) straightforwardly of all other quantities. Replacing every 



^ We remark that if there exists a velocity field v for which 9tB 4- V X (B X v) = 0, then 
magnetic flux is convected by v and the topology of magnetic field lines cannot change. In 
particular, if we merely add the Hall term to the ideal Ohm's law, then 9tB -|- V X (B X (u -|- 
m^—m^ J ^ B)), i.e., the magnetic field is essentially carried by the electrons. Hall-mediated fast 
reconnection requires a small amount of resistivity as well. 

To see this first take the divergence of Maxwell's evolutions equations. Then either (1) 
take the time derivative of Maxwell's constraint equations to eliminate the electromagnetic field 
and get a wave equation for the correction potentials, or (2) eliminate the correction potentials by 
taking the time derivative of the divergence of Maxwell's evolution equations and the Laplacian 
of the constraint equations to eliminate the correction potentials and get a wave equation for the 
divergence constraint error. 
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quantity X with nondimensional representation XXq gives a system governing the 
X values with exactly the same appearance, except that 1/e ~ c? . We drop hats. 

4. GEM magnetic reconnection challenge problem 

With the exception that our nondimcnsionalizcd light speed is 10 rather than 
their value of 20, our settings are equivalent to those of [J, which reflect the 
settings and conventions of the original GEM problem [3]. To map our Sl-like 
nondimensionalization onto their Gaussian-like nondimensionalization, rescale the 
electromagnetic field by Bgem = \/47rB and Eqem = -^^E. 

4.1. Computational domain. The computational domain is the rectangular 
domain [—Lx/2, Lx/I] x [— iy/2, ij,/2], where Lx = Stt and Ly ~ An. The problem 
is symmetric under reflection across either the horizontal or vertical axis. 

4.2. Boundary conditions. The domain is periodic in the x-axis. The bound- 
aries perpendicular to the y-axis arc thermally insulating conducting wall bound- 
aries. A conducting wall boundary is a solid wall boundary (with slip boundary 
conditions in the case of ideal plasma) for the fluid variables, and the electric field 
at the boundary has no component parallel to the boundary. We also assume that 
magnetic field runs parallel to and so does not penetrate the boundary (this follows 
from Ohm's law of ideal MHD, but we assume it holds gcnerally)Q 

4.3. Model Parameters. We carried out simulations for the following choices 
of the GEM model parameters: 



(1) 


rrii/me = 25, 


T,/T, 


= 5 


(original GEM), 


(2) 


mi/rrie = 1, 






and 


(3) 


mil me = 1, 




= 1. 





4.4. Initial conditions. The initial conditions are a perturbed Harris sheet 
equilibrium. The unperturbed equilibrium is given by 

'B{y) = Bo tanh(y/A)e:r, p{y) = 7r^n{y), 



2no 

ni{y) = ne{y) = no(l/5 + sech^{y/\)), Pe{y) = ^ p{y), 

X 

E = 0, p,{y) = jr-p^P{y)- 



On top of this the magnetic field is perturbed by 

(5B = — e^ X V(V'), where 
■4>{x, y) = -00 cos(27rx/La,) cos{Try/ Ly). 
In the GEM problem the initial condition constants are 

A = 0.5, Bo = l, m = l, 1^0 = So/10. 



^We remark that with the symmetries of the GEM problem at the conducting wall boundary 
may also be regarded as a symmetry conditions for a periodic boundary if the solution on the 
entire domain is reflected across its bottom boundary and negated, allowing the GEM problem to 
be solved on a doubled domain with periodic boundaries and infinitely smooth initial conditions. 
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5. Properties of the GEM problem 

5.1. Reconnected flux. In defining and discussing magnetic flux wc restrict 
ourselves to the first quadrant of the domain, referring to it as if it were the entire 
domain. 

The ideal MHD model implies the frozen-in flux condition, which says that 
magnetic flux is convectcd with the fluid. The boundary conditions and symme- 
tries of the problem dictate that at the boundaries fiuid can only move parallel to 
the boundaries and that the in-plane fluid velocity at the corners must be zero. 
Therefore, the frozen-in flux condition would say that the flux through any bound- 
ary must remain constant. 

In fact, for models which permit reconnection, the frozen-in flux condition 
breaks down near the X-point, allowing magnetic field to diffuse and allowing fleld 
lines to break and reconnect so that they pass through the horizontal axis. We 
therefore deflne magnetic reconnection to be the loss of magnetic flux through the 
vertical axis into the first quadrant: 

Definition 5.1. The reconnected flux Frccon is defined by 

FiMt) ■■= Bi dy, F,,,onit) := ^^ioft(O) - Fieft(i)- 

Jo 

Proposition 5.1. The rate of reconnection is minus the value of the out-of- 
plane component of the electric field at the origin (i.e. the X-point). 

Proof. 

dtF,,,on{t) = -dtFieft(i) = - r^^^ dtB^ dy = dyE^ dy = -E^iO), 

Jq Jo 

since E^ is zero at the conducting wall. □ 

5.2. Reflectional symmetries. The GEM problem has rcflectional symme- 
try across the horizontal and vertical axes. Wc impose this symmetry by restricting 
our computations to the first quadrant. At the origin, symmetries across the hor- 
izontal and vertical axes mean that vectors have only an out-of-plane component 
and pseudovectors must be zero. 

5.3. Reductions in Ohm's law at the origin. Symmetries at the origin 
reduce Ohm's law to: 

E3 =??J3 + -(me(9,,P.,l,3 + a,,P,,2,3) - mz(9.iPe,1.3 + 5:.,Pe.2,3)) 
P 

+ -^-^ ( at J3 -f J3V • u U3V • J -f — ^ J3 V • J 

p \ p 

If we neglect resistivity, off-diagonal pressure components, and assume that the fiow 
of charge and mass toward or away from the origin is small, then this becomes 

-dtF^ccon = E3 « ^'^^ dtJa ~ fh^rhedt{iz/ p), 
P 

i.e., the value of —ihirheSzl P at the origin will tend to track with i^rccon- 



^ This confirms the theoretical fact that an MHD model which only includes the B X u and 
Hall terms in Ohm's law cannot give fast reconnection, since both these terms must vanish at the 
origin. 

^ In simulations which incorporate the entire domain symmetry is often lost due to compu- 
tational noise and the inherent instability of the problem. 



TWO-FLUID MODEL OF ELECTRON-POSITRON PLASMA 



7 



6. Results 

Following the precedent of [§], [4], and [5], we carried out simulations of 
the GEM problem using a third-order shock-capturing Runge-Kutta Discontinu- 
ous Galerkin solver for a coUisionless two-fluid model with isotropic pressure for 
each species. (We enforced the divergence cleaning for the magnetic field but not 
for the electric field.) We carried out simulations on a quarter domain (hence 
enforcing symmetry) for mesh sizes of 32 x 16, 64 x 32, and 128 x 64. 

For the low mesh resolution we seemed to observe fast reconnection in the 
electron-positron plasma (based on the pattern of magnetic field lines and the rates 
of reconnection), but not for high resolutions. We hypothesize that the electron 
inertial term coupled with sufficient (numerical) resistivity is sufficient to yield fast 
reconnection. In future studies we hope to explicitly introduce collisional (diffusive) 
terms (rather than relying on numerical diffusion) and explore whether we can show 
convergence to fast reconnection for resistive isotropic two-fluid plasma. 

We plotted reconnected flux, --•fhifheJ^/ p., and a crude calculation of - Jj, 
(the cumulative sum of —E^ evaluated at integer times) for all three combina- 
tions of mass ratio and temperature ratio. Our plots indicate that —fhi'ftieJy,/ p 
only tracked with reconnected flux for the high-resolution simulations of electron- 
positron plasma up until the time when a large magnetic island formed around the 
origin. (These were the plots where fast reconnection did not occur.) 
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Reconnection for m./m^=1 , T./T^=1 (mx=128, my=64) 
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Figure 1. Three plots attempting to show convergence of reconnec- 
tion data for electron-positron plasma with equal initial temperatures. 
Reconnection is suppressed for finer mesh resolution but fast reconnec- 
tion occurs for coarse mesh resolution. We conjecture that numerical 
difi^usion coupled with the ion inertial term is sufficient to admit fast re- 
connection. The appearance of a large central magnetic island beginning 
around time 15 is curious. 
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Reconnection for m./m^=1 , T./T^=5 (mx=128, my=64) 
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Figure 2. Convergence study of reconnection for electron-positron 
plasma with initial temperature ratio of 5. The reconnection rate would 
appear to decrease with decreasing numerical diffusion. Only for the 
finest mesh does a central magnetic island form. 
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Reconnection for m./m =25, T./T =5 (mx=32, my=16) 
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Figure 3. Convergence study of reconnection for electron-positron 
plasma with mass ratio of 25 and initial temperature ratio of 5. Fast re- 
connection occurs for both coarse and fine mesh resolutions. The anoma- 
lously early cessation of reconnection in the simulation at intermediate 
resolution is evidently related to the appearance of a large central mag- 
netic island beginning around time 20. 



